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The ADM formalism for two-point-mass systems in d space dimensions is sketched. 
It is pointed out that the regularization ambiguities of the 3rd post-Newtonian ADM 
Hamiltonian considered directly in ci = 3 space dimensions can be cured by dimensional 
continuation (to complex d's), which leads to a finite and unique Hamiltonian as d — > 3. 
Some so far unpublished details of the dimensional-continuation computation of the 3rd 
post-Newtonian two-point-mass ADM Hamiltonian are presented. 
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1. Introduction 

The problem of finding the equations of motion (EOM) of a two-body system within 
the post-Newtonian (PN) approximation of general relativity is solved up to the 
3.5PN order of approximation for the case of compact and nonrotating bodies [by 
nPN approximation we mean corrections of order (w/c)^" ~ (Gm / {rc^))" to New- 
tonian gravity] . The 3PN level of accuracy was achieved only recently. There exist 
two independent derivations of the 3PN EOM using distributional (Dirac delta's) 
sources: either ADM-Hamiltonian-based,^'^ or harmonic-coordinate-based. There 
also exists a third independent derivation of the 3PN EOM in harmonic coordinates 
using a surface-integral approach.^ 

To cure the self-field divergencies of point particles it is necessary to use some 
regularization method. It turned out that different such methods applied in d = 3 
space dimensions were not able to give unique EOM at the 3PN order. Only by em- 
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ploying dimensional continuation was it possible to obtain unambiguous results.^''* 
In this note we review the dimensional-continuation-based derivation of the 3PN 
two-point-mass ADM Hamiltonian. 

2. ADM formalism for 2-point-mass systems in d space dimensions 

We use units such that c=167rG'(j+i = l. Wc work in an asymptotically flat (d + 1)- 
dimensional spacetime with Minkowskian coordinates x=(x^, . . . , x'^). Particles 
are labeled by the index a G {1, 2}; masses, positions, and momenta of the particles 
are denoted by rua, Xa=(a;^, . . • ,a;^), and Pa={pai, • • • ,Pad), respectively. We also 
define: Ta := x - Xa, Va := |ra|, n^ := ra/r„; ri2 := Xi - X2, ri2 := |ri2| (|v| 
means here the Euclidean length of the d- vector v). The canonical variables of 
the theory consist of matter variables {xa,Pa) and field variables (7y,7r*^), where 
the space metric jij is induced by the full space-time metric on the hypersurface 
a;'^=const; its conjugate tt*-' can be expressed in terms of the extrinsic curvature of 
that hypersurface. 

Source terms in the constraint equations written down for two-point-mass sys- 
tems are proportional to the d-dimcnsional Dirac delta functions (5(x — x,j). We use 
the ADM gauge defined by the conditions (TT = transverse-traceless) : 



The field momentum tt*-' splits into a TT part tt^^. and a rest tt^^ (traceless but 
expressible in terms of a vector), tt'^ = n''^ + Tr^rj,. If both the constraint equations 
and the gauge conditions are satisfied, the ADM Hamiltonian can be put into its 
reduced form: 



The FN expansion of the reduced Hamiltonian is worked out up to the 3.5PN order: 



H = Y,ma + H^ + ifiPN + H2PN + H2.5FN + ^^SPN + i^s.sPN + 0{{v/cf). (3) 



3. Dimensional regularization of the 3PN Hamiltonian 

In Rcfs. 1 it was shown that the Riesz-implemented Hadamard regularization of 
the 3PN two-point-mass Hamiltonian performed in c? = 3 space dimensions gives 
ambiguous results. The ambiguities were parametrized by two numerical coefficients 
called ambiguity parameters and denoted by Wkinetic and ^static 

Dimensional continuation consists in obtaining the 3-dimensional Hamiltonian 
as lim^^a if3PN(rf), where H3P^{d) is the Hamiltonian computed in d space dimen- 
sions. This can be done straightforwardly if no poles proportional to 1/ (d — 3) arise 
when — > 3 (or if one shows that these poles can be renormalized away, as happens 
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in harmonic coordinates*). Reference 2 has shown that out of all terms building 
up the Hamiltonian density there are ten terms TA{d), A = 1^ . . . , 10, giving rise to 
poles when d ^ 3. It was checked that the poles produced by these terms cancel each 
other, thus lim^^s H3p^{d) exists. Moreover, it was shown that for all other terms 
the 3-dimensional regularization give the same results as dimensional continuation. 

Let be the 3PN Hamiltonian obtained in Refs. 1 by using an Hadamard 

"partie finie" (Pf) regularization defined in d = 3 space dimensions. To correct this 
Hamiltonian one needs to compute the difference AH^p^ :— lim^j^s _ff3PN(d) — 
^3PN- Only ten terms Ta contribute to AH^p^, therefore 

» 10 » 10 

Ai73PN = lim d^x J2 TA{d) - Pf d?x ^ T^(3). (4) 

A=l A=l 

Below we present three different methods which we used to compute AH^p^. The 
details of the 2nd and 3rd method were not published so far. Knowing AH^p^ one 
determines the values of both ambiguity parameters: Wkinetic = 41/24, Wstatic = 0. 

1st method. In Ref. 2 Ai/gpN was computed by means of the analysis of the 
local behaviour of the terms Ta around the particle positions x Xq . 

2nd method. It is possible to compute all c?-dimensional integrals in Eq. ^ 
explicitly. To do this one uses the Riesz formula 

.rf, ,P ./2 r((a + d)/2)T{{f3 + ci)/2)r( ~ {a + f5 + d)/2) 

12-7^ r(-a/2)r(-/3/2)r((a + /3 + 2d)/2) ' 

and the distributional differentiation of homogeneous functions, e.g., 

Pf((d-2) ^"-"- ''O ^ .wl .. M(x-x.). (6) 



dx^xi rt''^ \ ri ) dT{d/2~l) 

3rd method. Instead of d-dimensional Dirac distributions S one uses d- 
dimensional Riesz kernels 5^^ to model point particles: 

^(x-x„)=^lim<5.„(x-xj, '5e„(x-x.):=-I^|l|^r--. (7) 

Then one uses the formula ^ to calculate the integrals in Eq. (jH) and, at the end of 
the calculation, one takes the limit £i ^ 0, £2 — > 0. No distributional differentiation 
is needed. 

We have shown that these three methods yield the same final results. 
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